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Abstract 

Numerical evidence is presented for the existence of a new family of static, 
globally regular 'cosmological' solutions of the spherically symmetric Einstein- 
Yang-Mills-Higgs equations. These solutions are characterized by two natural 
numbers (m > l,n > 0), the number of nodes of the Yang-Mills and Higgs field 
respectively. The corresponding spacetimes are static with spatially compact 
sections with 3-sphere topology. 



There has been considerable progress in the study of the spherically symmet- 
ric Einstein- Yang-Mills (EYM) and Einstein- Yang-Mills-Higgs (EYMH) equa- 
tions, stimulated by the discovery of globally regular solutions of the EYM eqs. 
HI, for a review see e.g. |J. Up to now most of the attention has been fo- 
cused on asymptotically flat, smooth particle- like (self- gravitating sphalerons 
and monopoles) and black hole solutions. 

In this paper we present a new 'cosmologicaP type of globally regular so- 
lutions of the EYMH eqs. describing static, spherically symmetric spatially 
compact space-times. Analogous solutions ('static universes') have already 
been found in EYM theory || in the presence of a cosmological constant A, 
(EYMC). More precisely, numerical evidence has been presented in Ref. || in- 
dicating the existence of a discrete family of solutions indexed by the number m 
of nodes of the YM field, for a special set of values of the cosmological constant, 
{A(m), m = 1, . . .}. The m = 1 solution is particularly simple as its energy 
density is constant and it has a simple analytic form It appears that the 

presence of the cosmological constant is essential for the very existence of such 
solutions. Therefore at first sight it might appear surprising that such 'static 
universe' solutions also exist in an EYMH theory without a cosmological con- 
stant. In fact this should not come totally unexpected as the self-interaction 
potential of the scalar field, V(3>), can generate the necessary energy density 
to support a spatially compact space-time. Intuitively V(3>) can act as a 'dy- 
namical cosmological constant'. This viewpoint has been particularly stressed 
by Linde and Vilenkin || in their work on 'topological inflation'. 

The EYMH theory has three different mass scales, the Planck mass Mp\ = 
l/VG and the masses il% and Mh of the YM resp. Higgs field, giving rise to 
two dimensionless ratios a = M\y/Mpi and (3 = Mn/iV%. It turns out to be 
convenient to introduce an 'effective cosmological constant' A = a 4 /3 2 /4. Our 
(m = l,ra = 0,l,2,...) EYMH solutions bifurcate with the EYMC solution of 
|| for a = v3, • • • at A = 3/4 with a vanishing Higgs field. Varying a leads 
to 1-parameter families (n = 0, 1, . . .) with A determined as a (possibly multi- 
valued) function of a. Similarly for some discrete values a m ^ n our solutions 
bifurcate with the higher nodes cosmological EYMC solutions of Ref. || , when 
the Higgs field tends to zero. 

Following the notation of , we write the spherically symmetric line element 

as: 

ds 2 = e 2 ^dt 2 - e 2X ^dR 2 - r 2 (R)dQ 2 . (1) 
The 'minimal' spherically symmetric Ansatz for the YM field is 

W^T a dx^ = W{R){Txde + T 2 sm6dtp) + T 3 cos 9dy , (2) 

where T a denote the generators of SU (2) and for the Higgs field 

$ a = H(R)n a , (3) 

n a denoting the unit vector in the radial direction. The reduced EYMH action 
can be expressed as 



S = - / dRe {v+x ^ 



1 + e- 2A ((r') 2 + v'(r 2 )') - e' 2X r 2 Vi - V 2 , (4) 
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with 



and 



V 2 = - T^ + ~ {H 2 - a 2 ) 2 + W 2 H 2 . (6) 



2r 2 

Varying the reduced action (|j) one obtains the EYMH equations: 



1 _ e -2A U2 + v n r 2y 



1 + e~ ZA r' z - 2e~ A e~ A rr 



er x (r'e- X )' + e^"* f e^rv')' + < 



(r 2 e u ~ x H'y -e u+x H (^2W' 



) + 2e- 2X r 2 Vi - 2V 2 


= o, 


(7a) 


' -2e~ 2X r 2 Vi -2V 2 


= o, 


(7b) 


,-2x9^) +r dV 2 
dr dr 


= o, 


(7c) 




= 0, 


(7d) 


l + ^(H 2 a 2 )) 


= . 


(7e) 



Introducing the combinations 

N = e~ x r', k = e~ A (r' + rz/), (8) 
the field eqs. (^a-e) take the form: 

2kN = l + N 2 + 2W 2 + r 2 H 2 - 2V 2 , (9a) 
r = N, (9b) 
N W 2 

N=(k-N) 2 rH 2 , (9c) 

r r 



k = [i _ K 2 + _ ^ {H 2 _ a 2 )2 _ 2H 2 W 2 ]/r , (9d) 
r 2 



W = w(^ + H 2 )-( K -N^, (9e) 

H = 2H^ + ^(H 2 - a 2 )) - (« + JV)f . Of) 

where / := df/dcr = e~ A /'. Note that in Eqs. @ the remaining gauge freedom 
(i.e. diffeomorphisms in the variable R) is implicit in the choice of the inde- 
pendent variable a. A convenient choice (also for the numerical integration) is 
e A = 1. 

Next we recall that solutions with a regular origin at r = R have the 
expansion 

W[r) = l-br 2 + 0(r 4 ) , (10a) 
H{r) = ar + 0(r 3 ) , (10b) 

N(r) = 1 - (2b 2 + ^ + ^ r 2 + 0(r 4 ) , (10c) 
«(r) = 1 + [2b 2 - ^ - ^f- j r 2 + 0(r 4 ) , (lOd) 
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Figure 1: Solutions with m = 1,2 zeros of W and n = 0, 1 zeros of H. The 
figures display W, H (solid), r, k (dashed), and N (dotted) as functions of a. 



where a, b are free parameters. Solutions satisfying the conditions (10) which 
also stay globally regular contain (among others) asymptotically flat self- 
gravitating monopoles. Integrating the field equations (|9|) one finds that for 
the generic solution, however, N(a) becomes zero at some finite a = ctq with 
finite values of all the other dependent variables. As it is immediately seen from 
Eq. (§3) ./V(<7o) = implies stationarity of r at r(o"o) = rQ. Then in general 
r(a) decreases on some interval for a > ctq. We refer to such a point as an 
'equator'. In the case of the EYM eqs. it has been proven in Ref. |8| that for all 
solutions with an equator r(a) decreases from ro all the way to r = where a 
curvature singularity develops ('bag of gold'). Adding a cosmological constant 
radically changes this conclusion and there exists an infinite family of regular 
solutions with 3-sphere topology ||. 

Since for H(p) = the EYMH system reduces to an EYMC theory where 
the cosmological constant is given as A = a 4 /? 2 /4, it is natural to search for 
solutions of Eqs. (0) with a nontrivial Higgs field bifurcating with the EYMC 
ones. To see if such a bifurcating class really exists one should first establish 
the existence of regular solutions of the linearized Higgs-field equation (Qe) in 
the background of a 3-sphere type solution of the EYMC system since the 
bifurcation occurs for H — > 0. To discuss the linearization of the Higgs-field 
equation around solutions of the EYMC equations, it is more convenient to 



choose a different gauge from e 



p defined by dr/dp 



e x N. 



= 1, namely v = A and use the variable 
Then the linearization of the Higgs-field 



equation (^e), around an EYMC background solution (A = A&, r = rb(p), 
v = v b {p), W = W b {p), H = H b = 0) with A b = cfiffi/A) and H = h(p), reads 
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Table 1: Bifurcation points. 



m 


A m 


«m,0 


oc m ,x 


a m ,2 


a m ,3 


a m ,4 


1 


0.75 


1.732051 


0.707107 


0.480384 


0.369274 


0.301511 


2 


0.3642442 


1.929400 


1.080086 


0.692265 


0.521313 


0.420586 


3 


0.2932176 


2.694371 


1.331930 


0.957678 


0.712995 


0.569524 


4 


0.2703275 


3.629807 


1.686984 


1.170945 


0.922451 


0.739788 


5 


0.2608951 


4.730633 


2.061886 


1.424943 


1.102191 


0.912326 


10 


0.2512791 


12.441074 


4.080242 


2.777478 


2.137747 


1.745428 


20 


0.2501165 


37.982205 


8.256034 


5.578608 


4.284836 


3.494772 



as: 

(r!h')' = 2e^(wi-^h. (11) 
Let us consider first the simplest, m = 1 EYMC solution 0]: 

Wb = cosx, H b = , rf ) = V / 2sinx, (12a) 
Nb = n b = cos x , v h = , A b = 3/4 , (12b) 

where x = p/y/2. Then Eq. (|TT|) in the background (|l2|), reduces to the following 
simple equation: 

(sm 2 (x)h'(x))' = ^2cos 2 (a;) - -^sin 2 (x)^ h(x) , (13) 

where now ' stands for the derivative with respect to x. We note that Eq. (|l3|) 
can be transformed to a hypergeometric equation, however, the solutions of 
interest, i.e. regular on the background geometry and satisfying the condi- 
tion (|To|b) , can be directly found by the following trigonometric polynomial 
Ansatz for h{x) 

n 

h n (x) = Y,c k sm 2k+1 x. (14) 

k=0 

One then easily obtains the recursion relation for the coefficients c&: 
Ck (2fc-l)(2fc + l)-2-3/a 2 

Ck ~ (2k + 2)(2k + l)-2 Ck ~ U fe - 1 -" n > ^ 

together with the termination condition (c n+ i = 0) 

3 

a?. = - — r , n = 0, 1, . . . (16) 

1,n (2n + l)(2n + 3)-2' ' ' v ' 

yielding a\ = 3, a\ 1 = 3/13, a\ 2 = 1/H) • • • 

That is we have found an infinite family of regular solutions of the linearized 
Higgs field equation (bounded 'zero modes') in the m = 1 EYMC background. 
These 'zero modes', indexed by the number of zeros of h(x) (n = 0,1,...), 
indicate the existence of a family of globally regular solutions of the EYMH 
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Figure 2: The effective cosmological constant A = a 4 /3 2 /4 as a function of a 
for the solutions with m, n < 2. 

equations (with H{x) having the same number of zeros) bifurcating with the 
m = 1 EYMC solution for a = and a = a\^ n . 

From the o-dependence of the r.h.s. of Eq. (|i"3| ) it follows that the solutions 
h(x) vanishing at x = have n + 1 zeros for all a G [ai jra +i, ai )n )- According 
to standard theorems on Sturm-Liouville operators this number gives also the 
number of bound states of the differential operator in Eq. (|l3|). Thus whenever 
a crosses ai )n from above a new bound state appears, indicated by the existence 
of a zero mode for a = a\^ n . 

In contrast to the case m = 1, the higher node EYMC solutions {m > 2) 
are not known analytically. Nevertheless one can still conclude that for any of 
these EYMC backgrounds an infinite family of regular zero modes exists. To 
see this, it is sufficient to note that for a sufficiently small the r.h.s. of Eq. (|ll|) 
becomes arbitrarily negative over a finite interval of p. This implies that the 
number of bound states tends to infinity for a — > 0. According to the previous 
argument the same holds for the number of zero modes accumulating at a = 0. 
Numerical values for the first five bifurcation points with m = 1, . . . , 5, 10, 20 
are given in Tab. [l]. 

In order to verify the correctness of the above scenario we have integrated the 
field eqs. @ numerically. To simplify the singular boundary value problem, we 
have looked only for solutions, which are (anti)symmetric about the equator. 
This means we impose boundary conditions at r = and at the equator (a 
regular point). In addition to the vanishing of k at the equator, we require the 
vanishing of the functions W resp. W and H resp. H depending if m and n are 
even or odd. Thus 3 functions must have a common zero with ./V forcing us to 
tune 3 of the 4 available parameters a, (3, a, and b. Fig. |l| shows some solutions 
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with m = 1, 2 zeros of W and n = 0, 1 zeros of H. 

In Fig. ^ we have plotted the values of the parameters a and A for the 
solutions with m, n < 2. While a runs to large values with decreasing A for 
the n = solutions, i.e. with a nodeless Higgs field, the parameters for the 
solutions with n > 1 show the opposite behaviour. This structure remains true 
for higher values of m, although the change occurs in general for some uq > 1 
increasing with m (e.g. hq = 6 for m = 10). Similar to the case (2,0) some of 
the graphs in the a-A plane for m > 2 show a minimum of a, before they tend 
to large values of a (e.g. n = 1, 2, 4, 5 for m = 10). 

The families of numerical solutions for which a becomes large seem to have 
a smooth limit for a — > oo with A —* 1/4. In this limit, assuming that 
stays finite, Eqs. (|7|) reduce to an EYMH system with a cosmological constant 
A = 1/4 instead of a Higgs potential (/? = 0). The spatial sections of the 
corresponding space-times are no longer compact (a — > oo and ro — > \/2). We 
have also numerically integrated the limiting field eqs. and our results fully 
confirm the existence of the limit a — > oo. 

It was argued in [ O ] that the solutions of the EYMC system of pi] are 
unstable. As usual the criterion for instability is the existence of imaginary 
modes of the linearized time-dependent field equations. Although we do not 
quite approve of the methods of [11], we nevertheless believe that their result 
is correct. For the case m = 1 the instability has already been shown in ||. 
We expect the instability of the YM system to persist in our case with the 
Higgs field with the same number of unstable modes at least for small values 
of a. In view of the discussion of the solutions of Eq. (|Tl| ) given above there 
are additional instabilities of the EYMC solutions viewed as solutions of the 
EYMH theory (with a = 0) in the Higgs sector. More precisely, there are n + 1 
unstable modes for a G [ai m ,n+i> a m,n)- Turning to solutions with a 7^ but 
small the only questions is, if the number of the unstable Higgs modes of the 
(m,n) solutions is equal to n or n + 1. The answer depends on the behaviour 
of the zero mode at the bifurcation point a mj „ when a deviates from zero. If 
it turns into a bound state one gets n + 1, if it moves into the continuum n 
unstable modes. 

Similar to the globally regular solutions discussed in this paper, there are 
families of solutions with a horizon branching off from the corresponding so- 
lutions of |J without a Higgs field. In fact there are even more general so- 
lutions having two horizons, one of them replacing the regular origin of the 
forementioned class. We plan to give a detailed account of these solutions in a 
forthcoming publication. 
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